REPORT  DOCUMENTATION  PAGE 

'  form  ApQTOvrd 

0MB  No.  0704.0188 

»*vOfK  *<X  ol  *s  »0  —  »— t-,.  to'  ir.Mfs.n.o«x.  0«<4  vjwto. 

Ihf  datJ  j»'<*  tomcHrHfwj  .^r-d  |M^  ©«’  amr  avorn  o»  tK.i 

toi»--n>o«  ol  .rtOiaO«x|  »0^  ‘rCvKieKi  |h.^  Ov/dew  |0  (cf  ia\f<>.matiOrt  Oc>C<J«n>*vs  ♦aKj  >-oo-«v  »  /  U  JeMt-rvsrt 

n*»»»  1/04.  AiUnaitOA.  VA  jr-a  lO  ir'<'  «<  M  4  •«*.»<  IvHS-Kt  f  {0 /fl4  ^  H  »).  VVa^A-MjJOA.  OC  rtlO  J. 

1.  AGENCY  USE  ONLY  (lEJvf  0/jnA) 

2.  REPORI  DATE 

15  August  95 

3.  R£PORr  TYPS  ANO  DATES  COVtRfO  ! 

Final  Report  15  AUG  1994  -  15  AUG  1995 

4.  TITL£.AN0-SUBTITL£  '  ‘  -  ‘ 

3.  fUNOlMG  NUMBERS 

Dynamics  of  Cracked  Composite  Material  Structures 


5.  AUTHORIS) 

W.M.  Ostachowicz,  M.  Krawczuk,  A.  I^ak 

^  dtic 

^ELECTi: 

7.  PERFORMING  ORGANIEATION  .NAME(S)  AND  A00R£S5(ES) 

Polish  Academy  of  Sciences,  Institute  of  Fluid  Flow  Machinery 

14  Gen.  J.  Fiszera  Street,  80-952  Gdansk,  Poland 

8.  PERFOP 
REPORT 

F 

1 

9.  SPCNSCSING  .  MONITORING  AGENCY  .Nam£(5)  AND  AOOR£S5(ES) 

U.S.  Army  Research  Office 

Edison  House,  223  Old  Marylebone  Rd. 

London  NWl  5TH,  U.K 

:0.  SPC.MSCRING  .  .v.O.Ni;  O.Rl.MG  j 

AGE.MCY  .TEPC.RT  NUMSE.R  j 

N68171-94-C-9108 

n.  SUPPU'.'.U.NTARY  .motes 

17j.  OlST.TlOU  nCN  /  a  VAILAeiLiT  Y  S  T  A  T  E  .'.IE  .m  T 

Approved  for  pubhc  release;  distribution  unlimited. 

i:b.  oiS7Ri3u::c.w  czzi 

13.  ABST/^>4CT  2Q0\^or(ii) 


As  a  result  of  this  work  models  of  the  finite  beam  and  plate  elements  have  been  elaborated, 
to  enable  the  analysis  of  the  influence  of  the  fatigue  cracks  and  delaniinations  on  the  dynamic 
characteristics  of  the  constructions  made  of  unidirectional  composite  materials.  The  method  of 
modelling  the  crack  or  delamination  presented  in  the  report  enables  an  easy  modification  of  the 
elaborated  elements  according  to  its  specific  damage  (oblique  crack,  two-side  crack,  inside  crack, 
multiple  delaminations,  etc.). 

The  results  of  numerical  calculations  obtained  from  the  crack  model  are  in  consistence  with 
the  known  influence  of  the  position  and  depth  of  the  crack  on  the  decrease  of  the  natural  bending 
frequencies  of  the  cantilever  beam.  Simultaneously,  a  strong  influence  of  the  material  parameters  on 
these  changes  has  been  observed,  vdiich  does  not  exist  in  the  case  of  isotropic  materials. 

The  method  of  modelling  the  delamination  in  composite  beams  and  plates  is  versatile  and 
allows  analysis  of  the  influence  of  multiple  delaminations  on  natural  frequencies  of  beams  and  plates 
with  various  boundary  conditions.  Using  the  elaborated  models  effects  of  location  and  size  of 
delaraination  on  bending  natural  frequencies  of  composite  beams  and  plates  were  studied. 


14.  SUBUCT  TERMS 

Composite  Materials.  Craks  and  delamination 

15.  NUMBfR  OF  PAGtS 

t 

16.  PRICE  CODE 

\7.  SECURITY  CLASSIFICATION 

OF  REPORT 

UNCLASSIFIED 

IE.  SECURITY  CLASSiFICATIO.'l 

OF  THIS  PAGE 

UNCLASSFIED 

••9.  security  CLASSIFICATION 

OF  abstract 

UNCLASSIFIED 

20.  LIMITATIC.M  OF  ASSTRaC  . 

UL 

StjnCJfd  fotrr\  3*59) 

7r«^<r.O— <1 


19950911 


NS.N  75‘!O-01-jeO-S5OO 


Title  of  Research  Project : 

DYNAMICS  OF  CRACKED  COMPOSITE 
MATERIALS  STRUCTURES 


Name  of  Principal  Investigator : 

Prof.  Wiesfaw  Ostachowicz 


Name  of  Contractor : 

Mary  Ketelhut 


Contract  Number : 

N68171-94-C-9108 


The  report  number : 

5th  -  final  report 


Accesion  For 

1 

NTIS 

CRA&I 

S 

DTIC 

TAB 

□ 

Uiiann 

ounced 

□ 

Justification 

By 

Distribution/ 

Availability  Codes 

Dist 

Avail  and  1  or  I 

Spe 

cial 

The  report  period : 

15  August  1994  - 15  August  1995 


Ti\e  Research  reported  in  this  docum^J:^H^a^sdyeerun^e  possible 
th^gh  the  support  and  sponserf^p  of  the  U.C^'OovQrrment 
throi^  its  European  Reseapen  Office  of  the  U.S.  Army.  Thislt 
is  interfdedonly  for  the Jprternal  management  use  of  the  Contractor 
and  U.S^Gdv&mm 


Abstract: 

As  a  result  of  this  work  models  of  the  finite  beam  and  plate  elements 
have  been  elaborated,  to  enable  the  analysis  of  the  influence  of  the  fatigue 
cracks  and  delaminations  on  the  dynamic  characteristics  of  the  constructions 
made  of  unidirectional  composite  materials.  The  method  of  modelling  the 
crack  or  delamination  presented  in  the  report  enables  an  easy  modification  of 
the  elaborated  elements  according  to  its  specific  damage  (oblique  crack,  two- 
side  crack,  inside  crack,  multiple  delaminations,  etc.). 

The  results  of  numerical  calculations  obtained  from  the  crack  model 
are  in  consistence  with  the  known  influence  of  the  position  and  depth  of  the 
crack  on  the  decrease  of  the  natural  bending  frequencies  of  the  cantilever 
beam.  Simultaneously,  a  strong  influence  of  the  material  parameters  on  these 
changes  has  been  observed,  which  does  not  exist  in  the  case  of  isotropic 
materials. 

The  method  of  modelling  the  delamination  in  composite  beams  and 
plates  is  versatile  and  allows  analysis  of  the  influence  of  multiple 
delaminations  on  natural  frequencies  of  beams  and  plates  with  various 
boundary  conditions.  Using  the  elaborated  models  effects  of  location  and  size 
of  delamination  on  bending  natural  frequencies  of  composite  beams  and 
plates  were  studied. 
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1.  Introduction 

The  use  of  composite  materials  in  various  construction  elements  has  substantially 
increased  over  the  past  few  years.  These  materials  are  particularly  widely  used  in 
situations  where  a  large  strength-to- weight  ratio  is  required.  Composite  materials 
similarly  to  isotropic  materials  are  subjected  to  various  damages,  mostly  cracks  and 
delaminations.  They  result  in  local  changes  of  the  stiffness  of  the  element  made  of  such 
type  of  materials  and  consequently  of  its  dynamic  characteristics  are  altered.  Changes  of 
natural  frequencies  and  mode  shapes,  amplitudes  of  forced  vibrations  and  also  coupling 
of  vibrations  forms  are  observed.  The  dynamic  characteristics  of  damaged  elements  can 
be  correlated  with  the  location  and  size  of  damages  [1-4]  (cracks  or  delaminations). 
These  relations  are  frequently  used  in  diagnosis  of  such  constructional  elements  -  for 
example  [5-8]. 

The  problem  of  changes  in  dynamic  characteristic  of  constructional  elements  with 
damages  made  of  isotropic  materials  has  been  a  subject  of  many  papers,  the  review  of 
which  is  given  by  Wauer  [9],  but  only  a  limited  number  of  papres  have  been  devoted  to 
the  changes  in  the  dynamic  characteristics  of  composite  constructional  elements.  Adams 


1 


et  al.  [10],  found  that  damage  in  specimens  fabncated  from  fibre  reinforced  plastics 
could  be  detected  by  a  reduction  in  the  natural  frequencies  and  an  increase  in  damping. 
Cawley  and  Adams  [11],  successfully  tested  the  frequency  measurement  principle  on 
composite  matrix  for  unidirectional  composite  materials  in  the  presence  of  damage. 
Nikpour  and  Dimarogonas  [12],  presented  the  local  compliance  matrix  for  unidirectional 
composite  materials.  They  have  shown  that  the  interlocking  deflection  modes  are 
enhanced  as  a  function  of  the  degree  of  anisotropy  in  composites.  The  effect  of  cracks 
upon  the  buckling  of  an  edge-notched  column  for  isotropic  and  anisotropic  composites 
has  been  studied  by  Nikpour  [13],  He  indicated  that  the  instability  increases  with  the 
column  slenderness  and  the  crack  length.  In  addition  he  has  shown  that  the  material 
anisotropy  conspicuously  reduces  the  load-carrying  capacity  of  an  externally  cracked 
member.  Manivasagam  and  Chandrasekaran  [14]  have  presented  the  results  of 
experimental  investigations  upon  the  reduction  effect  of  the  fundamental  frequency  of 
layered  composite  materials  with  damage  in  the  form  of  cracks.  The  effects  of 
delamination  on  buckling  and  post-buckling  deformation  and  delamination  growth  with 
various  geometrical  parameters,  loading  conditions,  material  properties  and  boundary 
conditions  have  been  studied  extensively  by  Chai  et  al.  [15],  Bottega  and  Maewal  [16], 
Whitcomb  [17],  Yin  et  al  [18]  and  Chen  [19].  However,  only  a  few  investigations  have 
been  conducted  to  study  the  effect  of  delamination  on  vibration  characteristics.  Natural 
frequencies  of  delaminated  beams  have  been  studied  by  Ramkumar  et  al  [20]  on  the 
basis  of  the  Timoshenko  beam  theory.  The  authors,  however,  did  not  take  into  account 
the  effect  of  coupling  of  the  transverse  vibration  with  the  longitudinal  wave  motion  in  the 
upper  and  lower  split  layers.  Their  analytical  results  predicted  significant  reduction  of  the 
fundamental  frequency  (from  that  of  the  perfect  beam)  and  this  prediction  was  found  to 
disagree  with  the  experimental  observation.  Wang  et  al  [21]  used  the  classical  beam 
theory  but  they  considered  the  coupling  effect.  With  the  inclusion  of  coupling,  the 
calculated  fundamental  frequency  was  not  appreciably  reduced  by  the  presence  of  a 
relatively  short  delamination  and  the  results  were  in  close  agreement  with  experimental 
measurements.  The  effect  of  delamination  upon  the  natural  frequencies  and  mode  shapes 
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was  analysed  by  Shen  and  Grady  [22],  The  effect  of  coupling  between  longitudinal  and 
bending  vibrations  was  considered  in  their  model.  Free  vibration  of  a  laminated 
composite  beam-plate  with  a  one  dimensional  delamination  with  respect  to  postbuckled 
reference  states  was  also  studied  by  Yin  and  Jane  [23], 

It  is  characteristic  for  all  cited  papers  that  their  authors  applied  a  continues  or 
discrete  continues  models  i.e.  models  which  are  based  on  differential  equations.  The 
crack  in  such  models  was  substituted  by  a  spring  with  additional  boundary  conditions 
[12-13]  or  by  reduction  in  elastic  modulus  of  material  -  see  for  more  details  [10-11].  The 
delamination  was  modelled  by  additional  boundary  conditions  -  see  for  more  details  [22]. 
This  method  of  modelling  is  inconvenient  for  more  realistic  structures  than  beams  or 
plates  with  constant  cross-section.  For  this  reason  the  authors  of  the  presented  report 
elaborated  the  alternative  techniques  of  modelling  of  damaged  structures.  This  technique 
is  based  on  finite  element  method.  The  report  presented  the  formulation  of  characteristic 
matrices  for  beam  finite  element  with  the  transverse  crack  and  delamination  and  also  for 
the  plate  finite  element  with  delamination.  The  results  of  numerical  calculations  received 
on  the  basis  of  elaborated  models  are  also  presented  in  the  report. 


2.  Cracked  beam  finite  element 

2.1.  General  description  of  the  element 

In  the  Fig.2.1  a  composite  beam  finite  element  with  the  transverse  crack  is 
presented.  The  cracks  is  placed  in  the  middle  of  the  element  and  remains  open,  its  depth 
is  a.  The  width  of  the  element  is  B,  the  length  L  and  the  height  H.  The  angle  between  the 
fibre  and  the  axis  of  the  element  (plane  perpendicular  to  the  crack)  is  a.  The  element  has 
three  nodes  {I,  II,  III).  Each  of  them  has  two  degrees  of  freedom:  transverse 
displacements  q^,q^,qs  and  rotations  q2,q^,q6-  Considering  only  the  case  of  flat 
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bending,  and  assuming  that  there  is  no  warping  in  the  transverse  cross-section  of  the 
element,  displacements  on  the  both  sides  of  the  element  could  be  expressed  by  relations: 


Fig. 2.1  Composite  beam  finite  element  with  a  transverse  crack. 


\^^yi(x,y)  =  v,(x)  ’  \uy^{x,y)^v^{x) 

where  (/  =  1, 2)  denotes  rotation  and  V;  {i  -\,2)  -  transverse  displacement.  Transverse 
displacements  v,.  on  both  sides  of  the  crack  (on  left  and  right  hand  side  of  the  element) 
could  be  approximated  by  cubic  polynomials  while  the  independent  rotations  by 
quadratic  polynomials  [24]: 

I V,  (x)  =  ^-a^x  +  a^x^  +a^x^ 

[  (x)  =  a^+  a^x  +  UjX^ 

Assuming  that  the  distribution  of  the  shear  stress  is  linear  [25],  the  relations  (2.2)  could 
be  written  in  the  form: 

|vj(x)  =  «!  H-a^x  +  ajX^ -fa^x^  j^v^ix)  =  +a^x  +  a^x^ +a^QX^  (2  2  a) 

[{2^,  (x)  =  0:5  +a^x  +  3a^x^  “^11  +  3a,gX^ 

Using  the  conditions  in  the  nodes  of  the  element  and  the  conditions  expressing 

consistency  of  displacements  in  x  and  y  directions,  and  the  balance  of  forces  and 


I  2  3 

\v^{x')  =  +  UgX  +  a^^x  -t-a^x 

UJx)  =  a„+a^^xJra,,x‘ 


(2.2) 


moments  we  could  obtain  the  unknown  values  Values  and  represent 

flexibility  coefficients  of  the  element  in  the  crack. 


v,(x  =  0)  =  ^,  ,  ^i(^  =  0)  =  ^2 

V,  (x  =  1/2)  =  ^3  ,  (x  =  Z/2)  = 

Vj  {x  —  L^  —  q^  ,  (f>2  —  ^6 

~  -^Z^)  “  ^x2  ~  -^Z^)  —  '  ^x\  ~  -^/2) 

<,(x  =  L/2)  =  m;3(x  =  Z/2)  (2.3) 

=  Ll2)  =  Uy2i^  =  Ll2) 

w(,i(x  =  Ll2)  —  Uy2{x  =  Z-/2)  =  ^f^y\ ~  -^/2) 

<Xx=:L/2)  =  M"(x  =  i/2) 
w-(x  =  L/2)  =  M;'(x  =  i:/2) 


Applying  the  standard  finite  element  formula  the  coefficients  a,  -a,^  could  be  written  as; 


a,' 

■!  : 

>  =  A,< 

'^1' 

>  ^ 
> 

* 

’^1' 

1  > 

where  matrices  A,  and  A2  take  the  following  form: 
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a 
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aL 

2 
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4 
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2 

aL 

~j5L 

aL 
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aL 

0 

l{L-2c^) 

0 

4 

0 

L 

0 

0 
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0 

1 

3L-4c^ 

0 

4 

0 

0 

— 

U 

L 

where  a~  L- 2c^,  f3-L-c^. 

Next,  from  the  relations  between  the  coefficients  a,  -  a, 2  and  the  displacements  we  could 
easily  obtain  the  relations  for  the  shape  function  Nj  and  N2  in  both  parts  of  the  element; 


N,  : 


N,  : 


N,=X,-Aj,  N,=:X,-A, 
where  matrices  X,  and  X2  have  the  form: 


(2.7) 

(2.7.a) 


X,=X2 


0  0  0  -3x  y  -y  -xy 

1  X  x^  0  0 


The  strains  in  the  element  could  be  obtained  from  the  following  relations; 


^  ^  <7^1  (x) 

eA’‘.y)  =  -y^ 


Substituting  the  relations  (2. 2. a)  into  equations  (2.9)  and  using  the  relations  (2.4),  (2.5), 
and  (2.6),  the  strains  inside  the  elements  could  be  expressed  in  the  node-displacement 
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function.  It  means,  that  the  linear  strain-nodal  displacement  relation  matrices  B,  and 
have  the  form; 


1 

^7 

> 

.^12. 

Bi-Xi'Aj,  B2-X2-A2 

where  matrices  X;  and  Xj  denote: 

_  _  0  0  0  -6xv  0  -V 

X,  =  x,  = 

'  ^  0  1  2x  0  -1  -X 


(2.10) 


(2. 10.  a) 


(2.11) 


Finally,  the  matrix  of  inertia  and  stiffness  matrix  of  transverse  cracked  element  could  be 
written  in  the  following  form: 

M,  =  p\  N^'N.dJ^  +pj  = 

=  pAljx^,X,dV,A,+pAljxlX^dV^A^  =  (2.12) 

=  pAlMlA,+pAlMlA, 

K,  =  jB[DB,JF;+jB2^DB2c/F2  = 

=  A\\x\X,dV,A^+Al\xlX^dV^A^  =  (2.13) 

ft 

=  A[KlA,+A[KX 

where  matrix  D  is  the  stress-strain  relation  matrix  for  the  unidirectional  composite 
material.  Substituting  into  equations  (2.12)  and  (2.13)  relations  (2.5)  and  (2.6)  and  after 
multiplying  and  integrating,  the  inertia  and  stiffness  matrices  for  both  parts  of  the  element 
take  the  following  form: 
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0 


=  BHL 


0 
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Sym. 


0  0  0  0 
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2 
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0 

6 

4 

12 

0 

8 

16 

2 

8 

24 


(2.17) 

where  and  are  the  elements  of  the  stress-strain  relation  matrix  (see  Appendix  C). 


2.2.  The  model  and  the  algorithm  for  calculating  the  flexibility  in  the  crack 


The  flexibility  coefficients  of  the  element  due  to  the  appearance  of  the  crack 
could  be  obtained  from  the  Castigliano  theorem: 


.  _  ^'7/ 
“  dP,  dP^ 


(2.18) 


where  U  is  the  additional  elastic  strain  energy  of  the  element  caused  by  the  crack,  while 
P.  and  Pj  denote  the  independent  nodal  forces  of  the  finite  element.  The  additional 

elastic  strain  energy  in  the  case  of  cracks  existing  in  the  unidirectional  composite 
materials  [12]  equals  to; 


U  = 


\\dXkI+d,,±k„:^k,„ 

A  I  '=1  '=1  >=> 


+D,±Kl}^dA 


(2.19) 


where  A  is  the  surface  of  the  crack,  Kj^  (j  - 1,11  ,i- 1,2,...,  n)  are  the  stress  intensity 
factors  and  Z), ,  are  the  coefficients  depending  on  the  material  parameters. 
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A  =-0-5^22  Im[lA+lAi] 

A2 

D2  =  O.Sij,  Im[i'j  +52] 


(2.20) 


The  formulae  for  calculating  values  s'!,  and  are  shown  in  Appendix  B. 

In  the  general  case,  stress  intensity  factors  for  composite  materials  are  not  equal  to  the 
stress  intensity  factors  calculated  from  the  solution  of  the  crack  problem  of  the  same 
geometry  in  the  isotropic  material.  According  to  the  results  presented  in  the  paper  [26] 
these  factors  could  be  written  as: 

K,  =  a,4^F,(alH,iflLlH,c)  (2-21) 

where  cr.  denotes  the  stress  acting  in  the  crack,  a  is  the  depth  of  the  crack,  H  is  the 
height  of  the  element,  A  ihe  correction  factors,  which  consider  the  finite  dimensions 

of  the  element  and  properties  of  the  material  (anisotropy  of  the  material),  while  the 
material  parameters  A,  C,  are  given  in  paper  [26].  It  is  also  shown  there,  that  for 
^[xlIh  >  2  intensity  factors  for  transverse  cracks  in  composite  materials  could  be  given 

as: 


K,  =  a,4^F,(alH)Y^(C)  (2.22) 

where  is  the  correction  function  which  takes  into  consideration  the  anisotropy  of 

the  material  [26]. 

Finally,  we  obtain: 


(2.23. a) 

(2.23. b) 

where  from  [26]: 


A- 


0.752 +  2.02^+0.37(l-sin  Tjf 

Tj  COST] 


(2.24.  a) 
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(2.24,b) 


tan;/  0.923 +  0.1 99(1 -sin  77) 


rj  cost; 

i;  =  1  +  0.  l(  C- 1)  -  0. 0 1 6(  1)'  +  0 . 002(  4"- 1)' 


yjE\  i-E'2: 


:2L_y  -^22 

2G, 


(2.24.C) 

(2.24.d) 


and  ^  =  alH,r}=  najlH,  while  B  is  the  width  of  the  element.  The  method  of 
calculating  of  the  material  parameters  is  shown  in  Appendix  A. 


3.  Delaminated  beam  finite  element 


3.1.  General  description  of  the  element 


A  discrete  model  of  a  delaminated  part  of  the  beam  is  presented  in  Fig.3.1  The 
delaminated  region  is  modelled  by  three  beam  finite  elements  which  are  connected  at  the 
tip  of  the  delamination  by  additional  boundary  conditions. 


Fig.3.1.  The  delaminated  region  of  a  beam  modelled  by  finite  elements. 

The  layers  are  located  symmetrically  with  respect  to  the  x-z  plane.  Each  element 
has  three  nodes  at  x  =  -LI2,  x-Q,x-  LH.  At  each  node  there  are  three  degrees  of 
freedom  which  are  axial  displacement  q^,q^,gj,  transverse  displacement  q^,q6,qg  and 
the  independent  rotation  ■  Additionally  it  is  assumed  that  the  number  of  degrees 

of  freedom  is  independent  of  the  number  of  layers. 
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3.1.1.  Description  of  the  element  number/ 


Neglecting  warping,  the  displacements  u  and  v  of  a  point  can  be  expressed  as; 

\u{x,y)^u\x)-y-(t){x) 

lv(x,;/)  =  v°(x) 

where  m°(x)  denotes  the  axial  displacement,  ^(x)  the  independent  rotation,  and  v°(x) 
the  transverse  displacement. 

In  the  finite-element  modelling,  the  bending  displacements  v“(x)  are  assumed  to 
be  cubic  polynomials  in  x,  while  the  axial  displacement  m°(x)  and  the  rotation  ^(x)  are 
assumed  to  be  quadratic.  Additionally  it  is  assumed  that  shear  strain  variation  is  linear,  as 
proposed  by  Tessler  and  Dong  [25].  Employing  the  above  conditions,  the  displacements 
and  rotation  in  the  element  may  be  written  in  the  following  forms: 

M°(x)  =  aj  +a2X  +  a^x^ 

<  (l){x)  =  a^+a^x-{-'ia^x^  (3.2) 

v°(x)  =  ag  -f-a^x+a^x^  +agx'^ 

The  constants  a,  -a^  can  be  expressed  in  terms  of  the  element  degrees  of  freedom  by 
using  the  nodal  conditions  in  the  following  forms: 

w°(x  =  -LH)  - 
<j){x  -  -LI2)  -  q2 
v°  (x  =  -LH)  =  ^3 
u\x  =  Q)=q^ 

(f>{x  =  0)  =  q^  (3.3) 

v°(x  =  0)  =  gg 

w“(x  =  LIT)  =  q-j 
(j)  {x-  LIT)  =  q^ 
v\x^LIT)  =  q, 
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Finally  we  obtain; 


-^1+^7 

L 

2{q,-2q,+qi) 

Z" 


^4  =  ^6 


—6^2  q^^  ~  Zg-^ 

""""  6Z 

..  2iq^-2q,+q,) 

J2 

_  ^(^3  ~  ^^6  Zg) 

31" 


(3.4) 


Taking  into  account  Eqs.  (3.4)  and  Eq.  (3.2)  we  can  determine  the  matrix  of  the  shape 
function  for  the  single  layer  of  the  element. 


N  =  X-A 


(3.5) 


where  matrix  X  has  the  form: 


\  X  -y  -xy  0  0  0  -2x^y 

000  0  0  \  X  x^  x^ 


whereas  the  matrix  A  can  be  expressed  as: 
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0  0  0 
-00 


0  0 


“T  0  0 

0  0  0 

0  0- 

L 

0  0  0 

1  1 


0  0 


Employing  the  shape  function  matrix  for  the  single  layer,  we  can  determine  the  inertia 
matrix  of  the  whole  element  using  the  following  formula: 


j=R  J=R 


M.  =  Zm;  =  ZP/ jN’'NrfF,  =ZP;A’'  J X^XdV,  A 

j'=i  ;=i  Vj  >=>  yj 

=  2p;A^MfA 

;=i 


where  j  denotes  the  number  of  a  layer,  R  the  global  number  of  layers  in  the  element,  Vj 
the  volume  of  the  y-th  layer  of  material  and  Pj  the  density  of  the  z-th  layer. 

The  value  of  the  integral  in  Eq.  (3.8)  (for  the  y-th  layer)  can  be  expressed  in 
closed  form  as; 
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Jj 

a  0  —  a 

12 

—  a  0 
12 

r 

—  a 
80 


0  0 


0  0 


0 

^  A  • 


0  0 


0  0 


L 

-—P 

8 


Mf  =  BL\ 


■j  0  0 


a  0 


a  0 


- a  + - r 

448  80  J 


where  a  =  -H^,  -H],  /  =  -  H] . 

The  strains  of  the  single  layer  of  material  are  given  by  the  following  formulas: 


^  ^  ^u{x,y)  ^  du\x)  a(f){x) 

dx  dx  dx 

dii{x,y)  ,  ^(x,y)  ^\x)  ^ 


(3.10) 


Taking  into  account  relations  (3.2)  and  (3.4),  the  strains  in  the  single  layer  can  be 
expressed  as  a  function  of  nodal  degrees  of  freedom: 


^2 

B  . 


(3.11) 


where  matrix  B  equals 


B  =  X  A 


(3.12) 
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and  matrix  X  is  given  as; 


0  1  2x  0  0  0  0  -6xy 

0  0  0  -1  -X  0  1  2x  0 


(3.13) 


The  stiffness  matrix  of  the  whole  element  has  the  form; 


j=R  j^R 


K.  =  Ek;  =z  J =E A"! X^D.XrfF,  A : 


;=1  ;=1  y 


;'=i  K, 


=  Ea^k;a 


(3.14) 


where  denotes  the  matrix  which  describes  relations  between  stresses  and  strains  in  the 

7-th  layer  of  the  element  (see  Appendix  C). 

The  values  of  the  integral  in  Eq.  (3.14)  (for  the  j-th  layer  of  the  material)  can  be 
presented  in  closed  form  as; 


‘0  0  0 


0 


0  0 


S^^a  0  0 


c 

0  0 


o  r2  s 


K{=BL\ 


0  -5,, a 

12  6 
0  0 


— - — a  0 


S,,I}  „ 
t  2 

SnL^r  . 


(3.15) 


where  a  =  -H],  y  =  H] . 
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3.1.2.  Description  of  elements  number//  and  III 

In  order  to  connect  element  I  with  elements  11  and  III,  the  following  boundary 
conditions  are  applied  at  the  tip  of  the  delamination; 

w°(x)  =  m°(x)  +  Z3^z53(x) 

V°(X)  =  V°(X)  =  V3°(X) 

where  Z2  and  Z3  denote  distances  between  neutral  axes  of  elements  /-//  and  /-///, 
respectively  (see  Fig.3.1). 

Taking  into  account  relations  (3.16)  and  (3.2),  the  relationships  between 
constants  a,  -a^  for  the  above-mentioned  elements  can  be  evaluated  as: 

(3.17.a) 

af  -  a{  -  ,  af  =ci{- 

<  af  =  ^2  -  y2(^i  ,  ^^2  =  ^2  “  ■  1^-^) 

~  «3  - 3y2‘^9  >  ^3^  =a[- 3y^al 

where  the  superscripts  /,  II  and  III  denote  the  number  of  the  element  in  the  region  of 
delamination. 

The  shape  function  matrices  for  the  elements  number  II  and  III  will  have  the 
following  forms; 


N2  =X-A2 

(3.18) 

N3=X-A3 

(3,19) 

ai  =  af 


I  - 


^111 

d-j  —  Cl~j  —  Cl-, 


ai  -  ai^  =  . 


al=ag  ^  af 
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where  A,,  (i  =  1, 2) 


A,.= 


0 


]_ 
~  L 
2^ 

1} 

0 

0 

0 

0 


0 

0 


0 

0 

0 

0 

0 

0 

\_ 

~  L 
2^ 

Jj 

0 


0 

fi. 

L 

2^ 

Jj 

0 

~  L 
0 

6 

0 

2 


1 

0 

1} 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

1 

0 

1} 

0 


0 

IfL 

J} 

1 

0 

0 

]_ 

3 

0 

4 


0 

L 

2^ 

1} 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

J_ 

L 

2^ 

1} 

0 


0 

_iL 

L 

0 

J_ 

L 

0 

_j. 

6 

0 

2 

3^. 


(3.20) 


Taking  into  account  matrices  A2  and  Aj  we  can  determine  (using  relation  3.8)  the 
inertia  matrix  of  elements  //  and  III. 

In  a  similar  way  matrices  B2  and  Bj  of  elements  II  and  III  can  be  evaluated,  and 
finally  the  stiffness  matrices  (using  relation  3.14)  of  these  elements  can  be  calculated. 


4.  Delaminated  plate  finite  element 

4.1.  General  description 

The  way  of  modelling  delaminated  region  in  a  composite  plate  with  delamination 
is  presented  in  Fig.4.1.  The  delamination  is  modelled  by  three  plate  finite  elements.  In 
order  to  connect  these  elements  in  the  delamination  crack  tip,  additional  boundary 
conditions  are  applied.  Material  layers  in  an  element  are  located  symmetrically  with 
respect  to  x-y  plane.  Each  element  has  eight  nodes  with  five  degrees  of  freedom. 
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Fig.  4. 1 .  The  region  of  delamination  in  a  plate  modelled  by  finite  elements. 


Axial  displacements  u,  v  and  w  in  a  single  layer  can  be  expressed  as: 

u{x,  y,z)^u\x,y)-z- (j)S^,y) 

<  v[x,y,z)  =  v\x,y)-z-<l)y[x,y)  (4.1) 

w(x,>/,z)  =  w°(x,y) 

where  denote  mid-plane  displacements,  while and  (j)y  denote  independent 

rotations.  To  approximate  the  axial  mid-plane  displacements  and  rotations  biquadratic 
shape  functions  for  eight-node  element  have  been  used. 


=  N3,,40-^ 


1^40  J 


(4.2) 
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iV,fe  t;)  = -1/4  +  f’ - 17;+ ir +2f- 17-2#^ 

N^(l4=V2  +  ri-l^-4^ri 

N,{^,  tj)  = -1/4 +  ^^+4ri+  if  +24‘T,+247f 

N,(4,ri)  =  V2-4-2rf+44rf 

N,(4,ji}=V2  +  4-2tf~44rf  (4.2,a) 

>/)  =  - V4  +  +  4r,+  if  - 2/fr]-24if 

N,(4n)  =  y2-v-2^+^(n 

N,[4  i?)  =  - V4  +  - 4rj+  if  -2^ij+24rf 


where  ^  =  xlL  and  77  =  ylB . 

Using  standard  finite  element  formulae  the  inertia  matrix  of  the  whole  element  can  be 
determined: 


M,  =  'f  M;  =  (4,3) 

;=1  ;=1  Vj 

where  j  denotes  the  number  of  the  layer,  R  total  number  of  layers  in  the  element,  F,  the 
volume  of  the y-th  layer  of  material  and  Pj  the  density  of y-th  layer. 

The  strains  in  the  single  layer  of  material  can  be  calculated  from  the  following  relations; 


0  du 

g  -  g  =— - 

"  *  "  dx  dx 

0  ,  ^‘t>y 

=  s^  +  z-  = - z - 

^  dy  dy 


r. 


^y  dx  J 


<^y  dx 


7  yz 


dw^ 

dx 

dw^ 

dy 


(4.4) 


Taking  into  account  relations  (4.1),  (4.2)  and  (4.4)  the  strain  in  a  single  layer  can  be 
expressed  in  the  following  form: 


20 


Thus  the  stiffness  matrix  of  the  whole  element  can  be  written  as: 


(43) 


K.='fK;='f  (43) 

;=i  ;-i  Vj 

where  D^.  denotes  the  stress-strain  relations  matrix  for  the  >th  layer  of  material  (see 
Appendix  C). 


4.2  Boundary  conditions  in  the  delamination  crack  tip 


To  connect  the  elements  /  with  elements  II  and  III,  modelling  the  delamination  area  and 
to  satisfy  continuity  of  the  displacements  the  following  conditions  must  be  fulfilled  (see 
Fig.4.2): 


Fig. 4. 2.  The  cross-section  of  the  plate  in  delamination  crack  tip. 


vl^vl-zl.6,  vl=v^,-zl.<l>^ 


where  2°  and  z®  are  the  distances  between  neutral  axes  of  elements  /-//  and  I-III. 
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Applying  relations  (4.7)  to  (4.1),  the  displacement  fields  of  elements  II  and  III  which 
model  delamination  region  are  obtained.  Similarly,  the  inertia  and  stiffness  matrices  of 
these  elements  can  be  calculated. 


5.  Numerical  calculations 

5.1  Natural  frequencies  of  cracked  composite  beam 

The  examples  of  numerical  calculations  showing  the  influence  of  the  crack 
parameters  (depth  and  position)  and  material  parameters  (volume  of  the  fibre  and  the 
angle  of  the  fibre)  on  the  changes  of  the  frequency  of  natural  bending  vibrations  were 
carried  on  for  the  cantilever  beam,  which  geometrical  dimensions  are  shown  in  Fig.  5.1. 
It  has  been  assumed  that  the  beam  is  made  of  unidirectional  composite  material  (graphite 
-fibre  reinforced  polyamide).  Materials  parameters  of  the  components  and  relations  to 
calculate  the  gross  material  coefficients  for  the  analysed  composite  material  are  presented 
in  Appendix  A. 


Fig.  5.1.  Geometry  of  the  cantilever  composite  beam  with  a  transverse  crack. 


22 


First  natural  frequency  [1/s] 


In  all  cases  we  assumed  ten  beam  finite  elements  presented  in  the  paper 


(including  the  element  with  the  crack)  to  divide  the  analysed  beam. 

The  numerical  calculations  were  carried  on  a  PC  computer  using  a  program 
written  by  the  authors.  For  the  standard  eigenvalue  problem  the  QL  method  was  applied. 


Fig.  5.2  Influence  of  the  crack  depth  on  the  first  three  natural  frequencies  of  the  cantilever 
composite  beam  a)  first  frequency,  b)  second  frequency,  c)  third  frequency. 


23 


First  natural  frequency  [Ms] 


Fig.  5. 3  Influence  of  the  crack  location  on  the  first  three  natural  frequencies  of  the 
cantilever  composite  beam  a)  first  frequency,  b)  second  frequency,  c)  third  frequency. 


Angle  of  fibres  [deg] 


Fig.  5.4  Influence  of  the  fibres  angle  on  the  first  three  natural  frequencies  of  the 
cantilever  composite  beam  a)  first  frequency,  b)  second  frequency,  c)  third  frequency. 


Fig.  5. 5  Influence  of  the  fibres  volume  on  the  first  three  natural  frequencies  of  the 
cantilever  composite  beam  a)  first  frequency,  b)  second  frequency,  c)  third  fi'equency 


In  Fig.5.2  and  Fig.5.3  the  results  showing  the  influence  of  the  depth  and  the 
position  of  the  crack  on  the  first  three  natural  frequencies  of  the  analysed  beam  are 
shown.  Numerical  calculations  have  been  carried  on  assuming  the  overall  volume  of 
fibres  of  20%  and  the  angle  of  the  fibre  45  degrees  (measured  between  the  geometric 
axis  of  the  beam  and  the  material  principal  axes). 

From  the  Fig.5.2  and  5.3  it  is  clear,  that  the  increase  of  the  depth  of  the  crack 
causes  in  each  case  decrease  of  the  each  calculated  natural  bending  frequency  of  the 
analysed  beam.  The  influence  of  the  position  of  the  crack  on  the  changes  of  the  /-th 
natural  bending  frequency  of  the  analysed  beam  is  more  complicated  and  should  be 
considered  together  with  the  character  of  the  vibrations  corresponding  to  that  frequency. 
It  shows,  that  the  decrease  of  the  considered  natural  frequency  is  most  substantial  when 
the  crack  is  placed  in  the  node  and  least  substantial  when  the  crack  is  placed  in  loop  of 
vibrations.  This  can  be  explained  by  the  fact  that  in  the  case  of  vibrations  of  the  cantilever 
beam  maximal  bending  moment  appear  in  the  nodes,  while  in  the  loops  of  vibrations 
bending  moment  is  assumed  zero.  The  decrease  of  stiffness  in  the  crack  corresponds 
(through  the  values  of  the  stress)  to  the  value  of  the  bending  moment.  Hence,  the  highest 
decrease  in  the  stiffness  caused  by  the  crack  (and  also  the  highest  decrease  of  the  natural 
bending  frequency  of  the  beam)  of  the  same  depth  will  occur  always  in  the  places  where 
the  bending  moment  is  highest. 

Fig. 5. 4  and  Fig. 5. 5  show  the  influence  of  the  material  parameters  (described  by 
the  angle  of  the  fibres  and  relative  volume  of  the  fibres)  on  the  changes  of  first  three 
natural  bending  frequencies  of  the  analysed  beam.  For  calculations,  we  assumed  non- 
changing  position  of  the  crack  {LJ L  =  Q.25)  and  constant  value  of  its  depth,  namely 

35%  of  the  height  of  the  beam. 

Two  conclusions  could  be  obtained  from  Fig.5.4  and  Fig. 5. 5.  Firstly,  with  the  increase  of 
the  angle  of  the  fibre  all  natural  bending  frequencies  of  the  beam  increase  with  respect  to 
the  values  of  the  natural  bending  frequencies  of  the  beam  without  the  crack.  For  all 
angles  greater  than  60  degrees,  the  decrease  in  the  natural  frequency  caused  by  the  crack 
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is  almost  unnoticeable.  This  conclusion  does  not  concern  the  case  of  isotropic  material, 
where  the  relative  volume  of  the  fibre  equals  1  or  0. 

Secondly,  as  a  result  of  the  numerical  calculations  it  has  been  shown,  that  the  decrease  of 
the  natural  bending  frequency  of  the  beam  caused  by  the  crack  strongly  depends  on  the 
volume  of  the  fibre.  These  facts  could  be  explained  by  noticing  the  changes  in  factors 
Dj ,  £>2  in  the  function  of  the  angle  and  relative  volume  of  the  fibre. 

5.2  Natural  frequencies  of  delaminated  beam 

The  formulation  of  the  elements  and  the  method  of  modelling  of  the  delaminated 
region  of  the  beam  have  been  evaluated  by  performing  several  example  calculations. 

Numerical  calculations  have  been  made  for  the  cantilever  beam  of  the  following 
dimensions:  length  600  mm,  height  25  mm  and  width  50  mm.  The  beam  was  made  of 
graphite-polyamide  composite.  It  was  assumed  that  all  layers  of  the  beam  have  the  same 
mechanical  properties,  i.e.,  the  volume  fraction  of  fibers  and  the  angle  of  fibers  in  each 
layer  are  identical.  The  mechanical  properties  of  the  applied  material  are  given  in 
Appendix  A. 

The  first  example  illustrates  the  influence  of  the  delamination  position  along  the 
beam  height  upon  the  changes  of  the  first  bending  natural  frequency.  The  length  of 
delamination  was  equal  to  120  mm  (a! L  =  0.2)  and  the  center  of  delamination  was 
located  300  mm  from  the  free  end  of  the  beam  (£,/£  =  0.5).  The  angle  of  fibers 
(measured  from  x-axis  of  the  beam)  was  45  degrees.,  whereas  the  volume  fraction  of 
fibers  was  equal  to  20%  the  volume  of  the  beam.  In  this  case  the  beam  was  discretized  by 
12  finite  elements  (4  elements  in  2  layers  around  the  region  of  delamination  and  8 
elements  outside  the  delaminated  region).  The  results  of  numerical  calculations  are  given 
in  Fig5.6.  It  is  clearly  shown  that  the  largest  drop  in  the  first  natural  frequency  is 
observed  when  the  delamination  is  located  along  the  neutral  axis  of  the  beam.  When  the 
delamination  is  located  near  the  upper  or  the  lower  surface  of  the  beam  the  changes  in 
the  first  natural  frequency  can  be  neglected. 
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37.0  ^ 


Relative  posKion  of  delamination  [h/H] 


Fig.  5. 6.  Influence  of  the  delamination  position  (along  the  beam  height)  on  the  first 
natural  frequency  of  the  cantilever  composite  beam. 

In  the  next  example  the  influence  of  the  length  of  the  delamination  upon  the  drop 
in  bending  natural  frequencies  of  the  analysed  beam  was  observed.  It  was  assumed  that 
the  delamination  growth  from  the  free  end  of  the  beam.  The  other  parameters  were  the 
same  as  in  the  first  example.  The  results  of  numerical  calculations  are  presented  in 
Fig. 5. 7.  It  is  noted  that  when  the  length  of  the  delamination  increases  the  values  of 
natural  frequencies  are  greatly  reduced.  The  intensity  of  these  changes  also  depends  on 
the  number  of  natural  frequencies  (i.e.,  the  mode  shape  and  the  location  of  delamination 
along  the  beam). 
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Fig.5.7  Influence  of  the  delamination  length  on  the  first  three  natural  frequencies  of  the 
cantilever  composite  beam  a)  first  frequency,  b)  second  frequency,  c)  third  frequency. 

The  third  example  shows  the  influence  of  the  location  of  delamination  along  the 
beam  on  the  drop  in  bending  natural  frequencies.  As  in  the  first  and  second  examples  the 
beam  was  made  of  polyamide-graphite  composite  material.  The  delamination  was  located 

along  the  neutral  axis  of  the  beam.  The  length  of  delamination  was  equal  to  37.5  mm 
(a/Z  =  0.0625).  Fig.5.8  illustrates  the  changes  of  analysed  frequencies  for  different 

locations  of  the  delamination.  It  is  clearly  shown  that  the  changes  in  natural  frequencies 
strongly  depend  on  the  location  of  delamination.  For  the  analysed  beam  the  largest  drop 
in  natural  frequency  is  observed  when  the  center  of  the  delamination  is  located  at  the 
node  of  mode  shape  associated  with  this  frequency. 
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Fig.  5. 8  Influence  of  the  delamination  location  on  the  first  three  natural  frequencies  of  the 
cantilever  composite  beam  a)  first  frequency,  b)  second  frequency,  c)  third  frequency. 

5.3  Natural  frequencies  of  delaminated  plate 

Numerical  calculations  for  the  cantilever  composite  plate  have  been  carried  out 
for  the  following  plate  dimensions;  length  240  mm,  width  120  mm  and  height  6  mm.  The 
plate  consisted  six  layers  of  materials  with  changing  angle  of  fibres  (+45  degrees  in  1,  3, 
5  layers  and  -45  degrees  in  2,  4,  6).  Each  layer  of  the  plate  was  made  of  graphite- 
polyamide  composite.  It  was  assumed  that  all  mechanical  properties  except  the  angle  of 
fibers  are  the  same  in  each  layer.  The  mechanical  properties  of  the  applied  material  are 
given  in  Appendix  A. 

The  first  example  illustrates  the  influence  of  the  delamination  position  along  the 
plate  height  on  the  changes  of  the  first  bending  natural  frequency.  The  length  of 
delamination  was  equal  to  30  mm  (a/Z.  =  0.125)  and  the  center  of  delamination  was 
located  105  mm  from  the  free  end  of  the  beam  {LJL  =  0.4375).  The  volume  fi'action  of 
fibers  was  equal  to  20%  the  volume  of  the  plate.  In  this  case  the  plate  was  discretized  by 
54  finite  elements  (42  elements  in  2  layers  around  the  region  of  delamination  and  12 
elements  outside  the  delaminated  region).  The  results  of  numerical  calculations  are  given 
in  Fig. 5. 9.  It  is  clearly  shown  that  the  largest  drop  in  the  first  natural  frequency  is 
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observed  when  the  delamination  is  located  along  the  neutral  plane  of  the  plate.  When  the 
delamination  is  located  near  the  upper  or  the  lower  surface  of  the  plate  the  changes  in  the 
first  natural  frequency  can  be  neglected. 
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Fig.  5. 9  Influence  of  the  delamination  position  (along  the  plate  height)  on  the  first  ^natural 
frequency  of  the  cantilever  composite  plate. 

In  the  next  example  the  influence  of  the  length  of  the  delamination  upon  the  drop 
in  bending  natural  frequencies  of  the  analysed  plate  was  observed.  It  was  assumed  that 
the  delamination  growth  from  the  free  end  of  the  plate.  The  other  parameters  were  the 
same  as  in  the  first  example.  The  results  of  numerical  calculations  are  presented  in 
Fig.5.10.  It  is  noted  that  when  the  length  of  the  delamination  increases  the  values  of 
natural  frequencies  are  greatly  reduced.  The  intensity  of  these  changes  also  depends  on 
the  number  of  natural  frequencies  (i.e.,  the  mode  shape  and  the  location  of  delamination 
along  the  plate). 
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Fig.  5. 10  Influence  of  the  delamination  length  on  the  first  three  natural  frequencies  of  the 
cantilever  composite  plate  a)  first  frequency,  b)  second  frequency,  c)  third  frequency. 


The  third  example  shows  the  influence  of  the  location  of  delamination  along  the 
plate  on  the  drop  in  bending  natural  frequencies.  As  in  the  first  and  second  examples  the 
plate  was  made  of  polyamide-graphite  composite  material.  The  delamination  was  located 
along  the  neutral  plane  of  the  plate.  The  length  of  delamination  was  equal  to  30  mm 
(uf/Z  =  0.125).  Fig. 5. 11  illustrates  the  changes  of  analysed  frequencies  for  different 

locations  of  the  delamination.  It  is  clearly  shown  that  the  changes  in  natural  frequencies 
strongly  depend  on  the  location  of  delamination.  For  the  analysed  plate  the  largest  drop 
in  natural  frequency  is  observed  when  the  center  of  the  delamination  is  located  at  the 
node  of  mode  shape  associated  with  this  frequency. 
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Fig. 5. 1 1  Influence  of  the  delamination  location  on  the  first  three  natural  frequencies  of 
the  cantilever  composite  plate  a)  first  fi-equency,  b)  second  frequency,  c)  third  frequency. 


5.  Conclusions 

As  a  result  of  this  work  models  of  the  finite  beam  and  plate  elements  have  been 
elaborated,  to  enable  the  analysis  of  the  influence  of  the  fatigue  cracks  and  delaminations 
on  the  dynamic  characteristics  of  the  constructions  made  of  unidirectional  composite 
materials.  The  method  of  modelling  the  crack  or  delamination  presented  in  the  paper 
enables  an  easy  modification  of  the  elaborated  elements  according  to  its  specific  damage 
(oblique  crack,  two-side  crack,  inside  crack,  multiple  delaminations  etc.). 
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The  results  of  numerical  calculations  obtained  from  the  crack  model  are  in 
consistence  with  the  known  influence  of  the  position  and  depth  of  the  crack  on  the 
decrease  of  the  natural  bending  frequencies  of  the  cantilever  beam.  Simultaneously,  a 
strong  influence  of  the  material  parameters  on  these  changes  has  been  observed,  which 
does  not  exist  in  the  case  of  isotropic  materials.  The  results  above  are  also  backed 
theoretically,  but  to  justify  them  it  is  necessary  to  carry  out  the  experimental  research  as 
well.  A  preparations  to  do  this  has  already  been  made. 

The  method  of  modelling  the  delamination  in  composite  beams  and  plates  is 
versatile  and  allows  analysis  of  the  influence  of  multiple  delaminations  on  natural 
frequencies  of  beams  and  plates  with  various  boundary  conditions.  Using  the  elaborated 
models  the  effects  of  location  and  size  of  delamination  on  bending  natural  frequencies  of 
composite  beams  and  plates  were  studied. 

Based  on  the  numerical  results,  the  following  conclusions  are  drawn: 

1)  The  delamination  in  cantilever  composite  beams  and  plates  causes,  as  expected, 
reduction  in  bending  natural  frequencies. 

2)  The  changes  in  natural  frequencies  are  a  function  of  the  location  and  length  of  the 
delamination. 

3)  When  the  center  of  delamination  is  located  at  a  point  where  the  bending  moment 
has  the  maximum  value  (for  analysed  mode  shape)  the  reduction  in  the  bending 
natural  frequency  associated  with  this  mode  is  largest. 

4)  The  largest  drop  in  the  bending  natural  frequencies  is  observed  when  the 
delamination  is  located  along  the  neutral  axis. 

5)  When  the  length  of  the  delamination  increases  the  drop  in  natural  frequencies  also 
increases. 
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Appendix  A 


The  properties  of  the  graphite-fibre  reinforced  polyamide  composite  analysed  in  the 
paper  are  assumed  as  follows  [9]: 


elastic  modulus 

A  =  2.756  GPa 

m 

=  275.6  GPa 

Poisson's  ratio 

v„  =  0.33 

Vj.  =  0.2 

rigidity  modulus 

G  =  1.036  GPa 

m 

=  114.8  GPa 

mass  density 

=  1600  kg/m^ 

Pj.  =  1900  kg/m^ 

I 

I 


The  material  is  assumed  orthotropic  with  respect  to  its  axes  of  symmetry  which 
lie  along  and  perpendicularly  to  the  direction  of  the  fibre.  The  gross  mechanical 
properties  of  the  composite  are  calculated  using  the  following  formulate: 


p  =  PfV+p„(\-v) 

En  --E^v+EJ^\-o) 

E,+E^-v[e.-E„)v 


F  =F  \  V  '  '  V^/ 

V|,  =  v/y+v„(l-u) 


V 


Gn  -  G„ 


1+  Vn,-ynE^lE,, 
1+ 

Gf  +  G^+[G^~G^v 
Gf+G^-[Gf-Gj)v 
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C  —  -^22 


where  u  denotes  the  volume  fraction  of  the  fibre.  The  principal  axes  1  and  2  are  in  the 
plane  of  the  composite  specimen  aligned  along  and  perpendicularly  to  the  fibre 
directions. 


Appendix  B 

The  complex  constants  5i  and  5^  are  the  roots  of  the  presented  characteristic  equation 

[12]; 


'2-b^^S  +(2&j2  +  2^1255  +  ^22  ^ 

where  constants  by  are  calculated  from  the  following  relations: 
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^11  —  2(6j2  +  ^55  ^  ^^22^ 

Ki  -  +2(^12  +^65)^'”^  ■^\^n' 

K^iK+K-  2^66  )w"«'  +  (m"  +  n" ) 

^6  =(^11-  bu  -  K  h'^  +  (*12  -  *22  +  ^6 
K  =  (*n  -  ^12  -  ^65)"'^  +(*12  -  *22  +  *66)"^'« 
*66  =  2(i„  -2Z),2  +*22  -  heWn^+K^im^  +«") 


where  m=  cos(a),  n-  sin(a),  {a  denotes  the  angle  between  the  fibre  direction  and  axis  of 
the  beam  perpendicular  to  the  crack-see  Fig.2.1).  The  terms  by  correspond  to  the 

situation  when  the  geometric  axes  of  the  beam  coincide  with  the  material  principal  axes. 
These  are  related  to  the  mechanical  constants  of  the  material  by: 


*n  = 


■^22 


fr 

l-Vn  — 

V  / 


+  (1-4) 


•"22 


*12  ~ 


■(l  +  '^23) 


*66  = 


Roots  of  the  characteristic  equation  are  either  complex  or  pure  imaginary  and  can  not  be 
real.  Thus,  the  four  roots  separate  into  two  sets  of  distinct  complex  conjugates.  The 
parameters  5,  and  correspond  to  those  with  positive  imaginary  parts. 


Appendix  C 


In  the  case  of  analysed  beam  elements,  the  stress-strain  relations  matrix  posses  the  form: 


■^.6  *^66 
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while  for  the  plate  element,  the  stress-strain  relation  matrix  can  be  written  as: 


Elements  of  matrix  D  can  be  calculated  from  following  relations: 

*^11  -  +^22^'^  +2(6’i2  +Sg^)m^n^ 

Su  =  (Su  +S22 -2SjmW  +n^) 

^22  -  +S^2^‘^  +2{S^2 

Si6={S^i-S,2-Sjm^n  +  {S,2-S22  +‘S’sJ«'w 

>^26  =  (‘^u  -  +  -*^22  +SeeWri 

See  =  2(^„  +S22-2S,2-See)mW  +See(m^  +«") 

S^=S^^m^  +*^55"^ 

Sa$  ~  (*^55  ~ 


where  m=  cos(a),  n=  sin(Q:),  {a  denotes  the  angle  between  the  fibre  direction  and  neutral 
axis). 

The  terms  corresponding  with  the  material  principal  axes  are  determined  in  all  cases 
by  following  formulas: 
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